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Question 1 [9 marks] 


(a) A study examined the lifetimes (hours) of 3 types of batteries (Type) operating under 


3 temperatures (Temp). 


Match each of the 3 interaction plots (1, 2 or 3) in Figure 1 to one of the following 


descriptions (i, ii, iii or iv). [3 marks] 


(i) Both Temp and Type effect but no interaction. 


) 
(ii) Interaction between Temp and Type. 
(iii) Type effect but no Temp effect. 

) 


(iv) Temp effect but no Type effect. 
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Figure 1: Interaction plots 


Question 1 is continued on page 3 
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Question 1 continued 


(b) 


Caterpillars go through “free growth” periods during each stage of their life. A bi- 
ologist is interested in checking to see if two different regression lines are needed 
to model the relationship between metabolic rate (Mrate) and body size (Bsize) of 
caterpillars for the growth periods and the no growth periods (indicator variable 


Fg = 1 for growth, 0 otherwise). 


State the model for predicting Mrate if the rate of change in the mean response with 
respect to body size is the same during the growth and no growth periods. Explain 


all terms. [3 marks] 


An experiment on wheat was laid out in a completely randomised design with three 
treatments (treat) and five plots per treatment. The variable of interest was yieldN, 
the crop yield (kg/plot). The yields from the same plots in the previous year (pre- 
treatment) were also recorded (X). The following model was fitted, with the table of 


estimates, standard errors and 95% Cls for the regression coefficients given below. 


Give a detailed interpretation of the estimate and the 95% CI given in the last line 


of the output. [3 marks] 
Table 1: 

mod<-1m(yieldN~treat*X) 
betaCI (mod) 
HHHHHHHHHHHE HH 

Estimate Std. Error 2.5 % 97.5 % 
(Intercept) -28.55 13.29 -58.62 1.53 
X 1.62 0.42 0.68 2.57 
treat2 20.57 15.97 -15.55 56.69 
treat3 18.71 16.56 -18.76 56.18 
treat2:X -0.62 0.53 -1.82 0.58 
treat3:X -0.38 0.55 -1.62 0.87 
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Question 2 [8 marks] 
Researchers collected samples of female lake trout from a lake in September and Novem- 
ber of three consecutive years, to investigate their fertility. One measure of viability of fish 
eggs is the percent dry mass (PctDM), which reflects the energy potential stored in the 
eggs. The age of the fish in years (Age) and the month (Month) in which the fish was 


observed were also recorded. 


(a) Explain in detail what the output in Table 2 implies. [2 marks] 


Cand.mod<-list() 
Cand.mod[[1]]<-lm(PctDM~Age*Month, data=FishEggs) 
Cand.mod[[2]]<-1lm(PctDM~Age+Month, data=FishEggs) 
Cand.mod[[3]]<-lm(PctDM“Age, data=FishEggs) 


Modnames<-c("modeli", "model2", "mode13") 


K AICc Delta_AICc AICcWt Cum. Wt LL 


model2 4 120 0.00 0.79 0.79 -55.1 

modeli 5 122 2.70 0.20 0.99 -55.1 

model3 3 129 9.28 0.01 1.00 -61.1 
Table 2: 


(b) Describe the comparative relationship of the lines for the preferred model. Refer to 
Table 3. [2 marks] 


(c) Using the estimates of the regression coefficients write down the equations of 


the lines. [2 marks] 
(d) Predict the PetDM for an 18 year old fish caught in November. [2 marks] 
Estimate Std. Error 2.5 % 97.5 % 
(Intercept) 39.519 0.7783 37.934 41.105 
Age -0.229 0.0629 -0.357 -0.101 
MonthSep -1.519 0.4234 -2.382 -0.657 
Table 3: 
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Question 3 [10 marks] 
The growth of spores was measured in 10 dilutions of a potato flour suspension. Five 
plates were placed in each of the 10 dilutions, and the number of positive plates showing 
spore growth was recorded. 
The response (proportion of positive plates) was analysed as a function of the natural 
logarithm of the dilution (Idose). 


(a) Does the model appear to fit the data? Justify your response with reference to 
Table 4. [2 marks] 


Analysis of Deviance Table 
Model: binomial, link: logit 
Df Deviance Resid. Df Resid. Dev P(>|Chil) 
NULL 9 41.212 
ldose 1 37.979 8 3.233 7.15e-10 


Table 4: Analysis of deviance Table for GLM 


(b) Using the R output below (Table 5) state and test the relevant hypothesis concerning 


the predictor variable. Give an informative conclusion. [2 marks] 


Coefficients: 

Estimate Std. Error z value Pr(>|zl) 
(Intercept) 4.1737 1.2522 3.333 0.000859 
ldose 1.6226 0.4571 3.550 0.000386 


Table 5: Summary Table for GLM 


(c) How many plates, out of a sample of 10, would be predicted to show positive spore 


growth, given a log dilution of -2.0? Show all calculations. [3 marks] 


(d) Find and interpret the LD75. Show all calculations. [3 marks] 
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Question 4 [17 marks] 
In a 1970’s study, terminal cancer patients were given ascorbate to see if it would 
prolong their lives. A designed experiment indicated that ascorbate may indeed prolong 
their lives. A second question of interest was whether the effect of ascorbate was the same 
for all patients regardless of the organ affected by cancer? The second study only focused 
on cancer patients receiving ascorbate. Five different organs were represented among the 
patients, with each patient only presenting one affected organ: breast, bronchus, colon, 
ovary or stomach. Researchers wanted to model the association between survival time and 


type of cancer. 


(a) An exploratory data analysis was performed. The summary statistics (mean and stan- 
dard deviation for survival time in days) are given in Table 6, and a boxplot of patient 
survival is given in Figure 2. What does this exploratory analysis suggest? 

[2 marks] 


Table 6: 
Breast Bronchus colon Ovary Stomach 
mean 1396 212 457 884 286 


Std. Dev. | 1239 210 427 = 1099 346 
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Figure 2: 


Question 4 is continued on page 7 
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(b) Which power transformation is suggested by the following boxcox plot (Figure 3)? 


Explain your response. 


[2 marks] 


log-Likelihood 
-210 -200 -190 -180 -170 -160 -150 -140 


Figure 3: 


(c) An appropriate transformation was applied and a linear model between the transformed 


survival time and the type of cancer was fitted. 


(i) With reference to Table 7 and Figure 4 (p. 8) comment on the validity of the 


model assumptions. 


[2 marks] 


Table 7: 
Breast Bronchus colon Ovary Stomach 
mean (transformed) 6.56 4.95 5.75 6.15 4.97 
Std. Dev. (transformed) | 1.65 0.95 0.99 1.26 1.25 


Question 4 is continued on page 8 


Residuals 
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Question 4 continued 
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(ii) The estimates and 95% Confidence Intervals for the regression coefficients are 


given in Table 8. YT is the transformed variable of Survival. What conclusions 


can you draw from the output? [2 marks] 
Table 8: 
mod2<-1m(YT~ Organ) 
Estimate Std. Error 2.5 % 97.5 % 
(Intercept) 6.559 0.360 5.84 7.280 
OrganBronchus -1.605 0.462 -2.53 -0.680 
OrganColon -0.809 0.462 -1.73 0.116 
OrganOvary -0.408 0.607 -1.62 0.806 
OrganStomach -1.591 0.490 -2.57 -0.611 


Question 4 is continued on page 9 


Question 4 continued 
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(d) Further questions of interest looked at comparing groups. 


(i) Derive the contrasts to address each of the following questions, labelling them C1, 


C2, C3 and C4 respectively. 


[4 marks] 


C1 What is the difference in mean survival for those patients with breast or 


ovarian cancer compared with the mean survival time for the other three 


groups of patients? 


C2 What is the difference in mean survival for those patients with breast cancer 


compared with those with ovarian cancer? 


C3 What is the difference in mean survival for those patients with bronchus 


cancer compared with those with colon and stomach cancer? 


C4 What is the difference in mean survival for those patients with cancer of the 


colon compared with those with stomach cancer? 


(ii) Verify that Cl and C2 are orthogonal, showing your calculations. [2 marks] 
(e) Interpret the following output (Table 9) in the context of the problem. [3 marks] 
Table 9: 


mod3<-1m(YT~C1+C2+C3+C4) 
> betaCI (mod3) 
Estimate Std. 


(Intercept) 5.676 
C1 0.226 
C2 0.204 
C3 - 0.135 
C4 0.391 


Error 
0.1607 
0.0701 
0.3033 
0.1213 
0.2202 


2:5. 4 9F 
5.3543 5 
0.0860 0 

-0.4029 0 
-0.3779 0 
-0.0499 0O 


5h 
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Question 5 [12 marks] 

The weights (kg) of poplar trees given different treatments were measured at different 
sites. The two sites differ in soil type. The four treatments were control (no treatment: 
N), fertilized (F), irrigated (I), and fertilized and irrigated (FI). 


Three models (mod.lme, mod2.lme and mod3.lm) were fitted: 


mod. 1me<-lme(wt~ treat ,random=~1|site/treat ,data=poplar) 


mod2.1me<-lme(wt~ treat ,random=~1|site,data=poplar) 


mod3.1m<-lm(wt~treat, data=poplar) 


(a) (i) Which of the two factors (treat or site) is being treated as a random 


effect? [1 marks] 
(ii) Explain why it is appropriate that the factor identified in (i) should be treated 


as random. [2 marks] 


(b) Interpret the interaction plot (Figure 5 and indicate what the plot suggests as a likely 


model. [2 marks] 
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Figure 5: 


Question 5 is continued on page 11 
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Question 5 continued 


(c) Examine the model comparison in Table 10 and interpret the output. 


[2 marks] 
Table 10: 
anova(mod.lme,mod2.1lme, mod3.1m) 
Model df AIC BIC loghik Test L.Ratio p-value 
mod.lme 1 7 74.09 85.18 -30.05 
mod2.lme 2 6 77.42 86.92 -32.71 1 vs 2 5.322 0.0211 
mod3.1m 3 5 76.06 83.98 -33.03 2 vs 3 0.642 0.4231 
(d) The (edited) R output for the final model is given in Table 11. 
(i) Using correct notation, derive the variance components. [3 marks] 


(ii) What proportion of the total variance is attributable to residuals (unexplained 


variability)? [2 marks] 


Table 11: 
Random effects: 


Formula: ~1 | site 
(Intercept) 
StdDev: 4.881e-05 


Formula: ~1 | treat %in% site 
(Intercept) Residual 
StdDev: 0.3813 0.4508 


Please remember: This examination question paper MUST BE HANDED IN. Failure 


to do so may result in the cancellation of all marks for this examination. Writing your 


name and number on the front will help us confirm that your paper has been returned. 
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